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We  study  the  collapse  of  ultracold  fermionic  gases  into  inhomogeneous  states  due  to  strong  dipolar 
interaction  in  both  2D  and  3D.  Depending  on  the  dimensionality,  we  find  that  two  different  types 
of  inhomogeneous  states  are  stabilized  once  the  dipole  moment  reaches  a  critical  value  d  >  dc'. 
the  stripe  phase  and  phase  separation  between  high  and  low  densities.  In  2D,  we  prove  that  the 
stripe  phase  is  always  favored  for  d  >  dc,  regardless  of  the  microscopic  details  of  the  system.  In 
3D,  the  one-loop  perturbative  calculation  suggests  that  the  same  type  of  instability  leads  to  phase 
separation.  Experimental  detection  and  finite-temperature  effects  are  discussed. 
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The  rapid  experimental  development  of  ultracold 
atomic  and  molecular  physics  has  opened  up  new  op¬ 
portunities  to  study  quantum  many-body  ^sterns  with 
electric  and  magnetic  dipolar  interactions  [ll-l9|.  An  im¬ 
portant  feature  of  dipolar  interaction  is  its  explicit  spa¬ 
tial  anisotropy  of  the  d^^-szS-type  when  dipole  moments 
are  aligned  by  external  fields.  For  the  fermionic  dipo¬ 
lar  systems,  “^^K-^^Rb  has  been  cooled  down  almost  to 
quantum-degeneracy  [ij.  Anisotropic  Fermi  liquid  theo¬ 
ries  of  the  single  particle  and  collective  properties  have 
been  investigated  [iMi.  Furthermore,  unconventional 
Cooper  pairing  structures  have  been  studied,  including 
the  Pz-channel  pairing  in  the  single  component  systems 
[17h22||  ,  and  the  competition  between  the  s  +  d-wave  sin¬ 
glet  and  the  Pz-wave  triplet  channels  [23l-[^.  In  par¬ 
ticular,  a  novel  pairing  state  of  the  s  +  ip  type  broken 
time-reversal  symmetry  has  been  pointed  out  |25l| .  More¬ 
over,  magnetic  dipolar  systems  of  fermions  have  also  been 
experimentally  realized  Exotic  states  of  the  ferro- 
nematic  Fermi  liquid  and  unconventional  magnetic  states 
have  been  predicted  [H,  ■ 

Fermionic  systems  can  spontaneously  break  transla¬ 
tional  symmetry  in  both  charge  and  spin  channels.  Many 
years  ago,  Overhauser  pointed  out  that  even  in  the  weak 
coupling  regime  of  the  interacting  electron  gas,  a  2k f 
spin-density  wave  state  always  wins  over  the  uniform 
paramagnetic  state  at  the  Hartree-Fock  level  [2^.  How¬ 
ever,  correlation  effects  may  suppress  this  instability  and 
such  a  state  has  not  been  experimentally  confirmed.  Re¬ 
cently,  quantum  liquid  crystal  phases  in  strongly  corre¬ 
lated  systems  have  been  intensively  studied,  particularly 
in  doped  Mott  insulators  [1^.  Stripe  ordering  has  been 
observed  in  high  cuprates,  other  transition  metal  ox¬ 
ides,  and  quantum  Hall  systems  at  high  Landau  levels 

BS. 

In  this  article,  we  study  the  instability  toward  the 
spontaneous  inhomogeneous  phase  in  dipolar  fermionic 
systems.  In  two  dimensional  (2D)  systems,  the  strongest 
density-channel  instability  occurs  at  non-zero  momen¬ 
tum,  which  drives  the  density  wave  states  under  strong 
dipolar  interactions.  This  effect  is  based  on  the  pecu¬ 


liar  feature  of  the  Fourier  transform  of  the  dipolar  in¬ 
teraction,  which  is  robust  against  microscopic  details. 
However  in  three  dimensions  (3D),  the  instability  at  the 
one-loop  level  occurs  at  zero  momentum,  thus  it  leads  to 
phase  separation  into  high  and  low  density  regions. 

We  consider  the  single-component  dipolar  Fermi  gas 
with  dipole  moment  aligned  by  an  external  electric  field. 
The  long-distance  physics  of  this  system  is  described  by 
the  following  Hamiltonian 

k  k,k',q 

with  being  the  energy  dispersion  relation  and  p  the 
chemical  potential.  In  2D,  the  dipolar  interaction  in  the 
momentum  space  V2d(^  takes  the  form  [T^. 

F2d(<?)  =  2Tr(fP2{cos9)[ - g]  -|-  nqef  sin^0cos2(^q,  (2) 

where  d  is  the  dipole  moment;  (j)q  is  the  azimuthal  angle 
of  the  momentum  g  and  9  is  the  angle  between  the  di¬ 
rection  of  the  dipoles  and  the  z-axis.  e  is  a  short-range 
cutoff  roughly  equal  to  the  thickness  of  the  system  along 
the  z-direction.  Equation  ([2])  is  valid  for  ge  <C  1.  How¬ 
ever,  as  shown  below,  the  main  conclusion  of  this  paper 
is  independent  of  the  value  of  V{q)  at  large  g.  Equa¬ 
tion  ([2|)  contains  an  isotropic  part  which  depends  on  the 
microscopic  cut-off  e  and  an  anisotropic  part  which  de¬ 
pends  on  the  azimuthal  angle  of  (fiq.  We  first  consider  the 
purely  anisotropic  case  at  9  —  9o  =  cos“^  and  then 
the  interaction  simplifies  into 

F2d(^  =  ^^gcos2(/),.  (3) 

The  major  feature  of  Eq.  ([3])  is  its  linear  dependence  on 
g,  which  play  important  roles  in  driving  the  instability  at 
non-zero  wavevectors. 

The  stability  of  a  system  against  density  fluctuations 
is  determined  by  the  static  susceptibility  of  density  fluc¬ 
tuations: 

n2D(g,w  =  0)  =  {p{q,UJ  =  0)p{-q,UJ  =  0)),  (4) 
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FIG.  1:  (online  color)  Diagrammatic  expansion  for  n2i^(^. 
The  solid  lines  are  free-fermion  propagators  and  the  dashed 
lines  represent  the  dipolar  interaction. 


where  p{q,uj)  =  Efc,a  is  the  density  oper¬ 
ator.  The  inverse  of  1120  =  0)  is  the  energy  gap 

^2d(^  of  creating  a  density  wave  at  momentum  k.  Us¬ 
ing  the  technique  of  the  diagrammatic  expansion,  A2d(^ 
can  be  computed  order  by  order  as 

A2d(^  =  1120(9,  w  =  0)  ^ 

=  [n™(g,a;  =  0)]-i  +  U2o(<7l,  (5) 

where  w)  is  the  one-particle-irreducible  density- 

density  correlation  [H,  [s^  shown  in  Fig.  [TJ  The  sta¬ 
bility  criteria  of  the  homogeneous  phase  corresponds  to 
the  energy  gap  A2o(<^  being  positive-definite  for  any  q. 
Otherwise,  the  creation  of  a  density  wave  could  lower  the 
total  energy  of  the  system  leading  to  the  condensation  of 
density  waves.  Obviously,  A2d(^)  is  positive-define  for 
any  q  in  the  noninteracting  Fermi  gas.  This  stability 
remains  for  a  weakly-interacting  Fermi  liquid.  On  the 
other  hand,  if  the  interaction  has  an  attractive  channel 
(e.g.  the  dipolar  interaction  in  Eq.  ([3])),  an  inhomoge¬ 
neous  state  could  take  advantage  of  the  attractive  inter¬ 
action.  Under  strong  enough  attractive  interactions,  i.e., 
V2D  —00,  the  homogeneous  ground  state  is  unstable. 

We  begin  from  the  weak  coupling  side,  by  tuning  up 
the  interaction  strength.  In  general,  there  is  a  critical 
dipole  moment  dc  such  that  A2d(<^  >  0  at  any  q  for 
d  <  dc,  but  A2d(^  <  0  at  some  9  for  d  >  dc-  At 
d  =  dc,  A2D(iQ)  =  0  for  the  momenta  ±<5  and  remain 
positive  at  other  momenta  [35l|.  If  we  increase  d  further 
above  dc,  the  density  wave  fluctuations  with  momentum 
q  ^  Q  undergo  an  instability  and  condense.  Depending 
on  whether  Q  is  0  or  not,  this  instability  has  two  differ¬ 
ent  fates.  If  Q  >  0,  the  condensation  carries  a  nonzero 
momentum,  which  leads  to  a  stripe  state:  a  unidirec¬ 
tional  charge-density- wave  state  with  wavevector  Q.  On 
the  other  hand,  if  Q  =  0,  the  condensation  takes  place 
at  an  infinite  wavelength  and  results  in  phase  separation 
between  high  and  low  density  regions.  Here,  the  typical 
size  of  these  high  (low)  density  regions  and  their  spatial 
arrangements  are  determined  by  the  microscopic  details 
of  the  system.  In  the  particular  case  of  a  dipolar  gas,  the 
phase  separation  is  sensitive  to  the  details  of  the  short- 
range  behaviors  of  the  interaction. 

Due  to  the  2D  space-inversion  symmetry  r  — >■  — r, 
0;)]“^  is  an  even  function  of  9.  As  a  result,  the 
Taylor  expansion  of  A2d(^  at  small  momentum  must 


FIG.  2:  (online  color)  The  stripe  phase  in  2D  dipolar  gases. 
The  three  dotted  arrows  mark  the  x,  y  and  ^  axes.  The  lighter 
(darker)  regions  in  the  2D  plane  represent  higher  (lower)  den¬ 
sity.  The  solid  line  with  an  arrow  represents  the  direction  of 
the  external  field,  which  aligns  the  direction  of  the  dipoles. 

take  the  following  form 

A2d(9)  =[niS'(9-'=  0)]-i  +  ^gcos2</., 

+  2  ^  [n2D^(^]  ^Ig=09j9j  +  ■  •  ■  (6) 

We  focus  on  the  most  unstable  direction  =  ^ ,  in  which 
the  attraction  in  V2d('^  is  the  strongest  and  the  linear 
q  term  in  Eq.  ([6])  has  a  negative  coefficient.  Hence  the 
point  9  =  0  is  a  local  saddle  point  of  A2d('^-  Therefore, 
as  d  increases,  A2d(^  at  some  finite  momentum  will  be¬ 
come  negative  before  A2d(9  =  0)  could  reach  zero.  Thus, 
Q  must  be  a  nonzero  value,  which  indicates  that  the  sys¬ 
tem  will  collapse  into  a  stripe  phase  for  d  >  dc,  instead  of 
phase  separating.  This  is  the  main  conclusion  of  this  pa¬ 
per.  It  is  noteworthy  that  this  conclusion  only  depends 
on  the  small-momentum  (i.e.  long-range)  behaviors  of 
the  interaction  V  (9).  Therefore,  the  stripe  instability  is  a 
universal  property  of  the  2D  dipolar  Fermi  gas  at  d  >  dc, 
insensitive  to  the  microscopic  details  of  the  system. 

More  importantly,  the  conclusion  of  Q  >  0  is  a  non- 
perturbative  result  which  remains  valid  to  any  order  in 
the  loop  expansion  of  Fig.  [TJ  The  only  assumption  re¬ 
quired  here  is  the  analyticity  of  n2®(9,  w  =  0)  as  a  func¬ 
tion  of  9  near  q  =  0.  In  a  Fermi  liquid,  it  is  well-known 
that  n^^^((f,  w  =  0)  is  non-analytic  at  9  =  2kF  where  kp 
is  the  Fermi  wavevector.  However,  there  is  no  reason  of  a 
non-analytic  behavior  at  small  momentum  in  any  Fermi 
liquid  to  our  knowledge.  The  higher-order  diagrams  in 
Fig.  [1]  contain  interaction  lines  which  are  non-analytic 
at  small  momenta.  However,  due  to  the  fact  that  dipo¬ 
lar  interactions  are  short-ranged  in  2D,  we  do  not  expect 
any  singularity  for  n2®((f,  w  =  0)  around  9  =  0. 

We  check  this  conclusion  using  the  loop  expansion  in 
Fig.  [Tjto  the  one-loop  level  for  a  continuum  2D  system 
=  k^ ! 2m  with  m  being  the  mass  of  the  particles)  and 
find 

A2d(^  =  — [ - , - h  ICOS20.,-],  (7) 

where  x  =  and  dc  is  the  critical  dipole  momentum 
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ory,  which  is  fundamentally  different  from  other  similar 
phases  at  d  0  [s^-  Due  to  this  unique  property,  the 
stripe  phase  at  0  =  0  is  unstable  at  any  finite  tempera¬ 
ture  due  to  thermal  fluctuations  [13,  ,  while  the  stripe 

phases  at  6*  0  has  a  power-law  quasi-long  range  or¬ 

der  blow  a  transition  temperature  Tkt  (See  below  for 
details). 

In  3D,  the  dipolar  interaction  is  V3d(^  = 

^^^P2{cos9^),  with  9q  measuring  the  angle  between  q 
and  the  dipole  moment.  Due  to  the  structure  of  V3D,  the 
linear  term  of  q  in  Eq.  ([6])  is  absent  in  3D,  and  hence 
we  could  not  exclude  the  possibility  of  Q  being  zero.  In 
fact,  the  one-loop  calculation  for  a  system  in  continuum 
shows  that 


,  .  1  ,  2m  |1  +  2;|s 

,,„(,T  =  _|(l  +  -(l-x^)log|  — I) 

+  ^(3cos^|-  1)], 

c 


(8) 


FIG.  3:  (online  color)  One-loop  calculation  of  A{q)  in  (a)  2D 
and  (b)  3D.  In  the  2D  case,  we  choose  ^  and  for  the 

3D,  9^  =  n/2.  From  top  to  bottom,  the  dipole  momentum  for 
each  of  the  curves  are:  d/dc  =  0,  0.9,  1,  and  1.1  respectively. 
The  dashed  vertical  line  marks  the  momentum  2fc_F. 


dc  =  y  2k^m  •  4'q  is  the  azimuthal  angle  of  the  momen¬ 
tum  q  and  ri{x)  is  a  step  function  with  ri{x)  =  1  for  a:  >  0 
and  r](x)  =  0  for  x  <  0.  As  shown  in  Fig.  |3(a)[  Q  —  2kF- 
If  higher  order  corrections  are  considered,  although  the 
value  of  Q  may  vary,  the  fact  that  Q  >  0  remains  valid. 

For  0  ^  9o,  the  same  conclusion  applies  for  the  single 
component  case.  The  e-term  in  Eq.  m  is  momentum  in¬ 
dependent,  which  corresponds  to  short-range  contact  in¬ 
teraction.  Thus  it  has  no  effect  on  the  single  component 
fermions,  and  V{q)  =  qnd^[—2P2{cos9)  -I- sin^  0  cos  2(()q] . 
Along  (/'q  =  f  which  remains  the  direction  of  the 
strongest  instability,  V (q)  is  negative  and  grows  linearly 
with  q  for  an  arbitrary  value  of  9,  and  the  above  conclu¬ 
sion  still  holds.  We  notice  that  even  for  9  =  0,  i.e.,  when 
the  dipoles  are  perpendicular  to  the  plane,  the  stripe 
instability  still  exists  at  large  values  of  d  although  the 
interaction  is  now  repulsive.  For  such  an  isotropic  sys¬ 
tem  (at  9  —  0),  it  is  often  assumed  that  three  charge 
density  waves  form  a  triangular  lattice  to  minimize  the 
breaking  of  the  rotational  symmetry  [s^.  However,  if 
the  rotational  symmetry  is  spontaneously  broken  before 
the  homogeneous  state  becomes  unstable,  a  stripe  phase 
can  be  stabilized  at  zero  temperature  [37j|.  In  this  case, 
the  stripe  ordering  spontaneously  breaks  both  the  rota¬ 
tional  symmetry  and  the  translational  symmetry  in  the 
direction  perpendicular  to  the  stripes,  while  only  the 
latter  is  spontaneously  breaking  in  the  stripe  phases  at 
9^0.  Due  to  its  unique  symmetry  breaking  pattern, 
the  low-energy  properties  of  the  stripe  phase  at  6*  =  0 
is  described  by  the  quantum  McMillan-de  Gennes  the- 


where  x  =  q/2kF  and  dc  =  {32kF'^^)-  As  can  be 

seen  from  Fig.  |3(b)[  <3  =  0,  implying  a  phenomenon  of 
phase  separation  for  d  >  dc-  As  discussed  above,  this 
is  sensitive  to  the  short-range  details  of  the  interaction 
and  no  universal  conclusion  is  available  for  this  phase 
separation. 

Now  we  briefly  discuss  the  multi-component  case  in  2D 
(e.g.  the  degrees  of  freedom  of  hyperfine  spins).  For  sim¬ 
plicity  we  only  consider  the  case  of  Al  =  2.  In  this  case, 
the  e-term  in  Eq.  ©  cannot  be  ignored,  which  naturally 
exists  in  the  inter-component  interactions.  Following  the 
same  argument  as  in  the  single-component  case  studied 
above,  the  instability  in  the  charge  channel  for  a  2D  sys¬ 
tem  with  9  >  9o  leads  to  a  charge  stripe  phase.  However, 
for  9  <  9o  the  short-range  part  becomes  repulsive,  so  the 
density  wave  state  is  not  favorable. 

The  stripe  phase  has  a  density  modulation  in  the  direc¬ 
tion  perpendicular  to  the  stripes,  which  can  be  detected 
directly  via  the  measurement  of  the  local  densities.  In  ad¬ 
dition,  any  scattering  experiments  would  show  two  inter¬ 
ference  Bragg  peaks  at  momentum  ±Q  in  a  stripe  phase 
due  to  this  density  modulation.  There  are  also  other  indi¬ 
rect  tools  to  detect  a  stripe  phase.  Considering  a  density 
wave  with  density  varying  along  the  y  axis  [e.g.  a  density 
profile  p{x,  y)  =  Pq  +  pi  cos{Qy  +  (p)  with  {x,  y)  being  the 
2D  coordinate  in  the  real  space],  if  we  place  this  system 
in  a  shallow  potential  trap  and  then  introduce  an  extra 
narrow  ID  potential  well  along  the  direction  of  the  stripe 
[e.g.  V{x,y)  =  it  is  energetically  favorable  to  have 
y  =  0  being  a  high  density  region  due  to  the  potential 
well  V  {x,  y).  If  we  move  the  potential  trap  along  the  y  di¬ 
rection  but  keep  the  location  of  the  potential  well  V  {x,  y) 
unchanged,  the  stripes  will  largely  remain  at  their  initial 
positions  if  the  displacement  of  the  trap  is  less  than  half 
the  wave-length  of  the  stripes  tt/Q,  and  the  stripes  will 
jump  over  a  distance  of  2tt jQ  if  the  displacement  be¬ 
comes  larger  than  tt/Q-  On  the  other  hand,  the  jump 
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would  not  happen  in  a  homogeneous  state.  This  jump 
will  lead  to  a  center  of  mass  oscillation  for  the  trapped 
particles  due  to  the  lack  of  dissipation,  which  could  serve 
as  an  indirect  signature  for  the  stripe  phase. 

Since  a  real  system  has  a  finite  temperature,  the  contri¬ 
bution  of  the  thermal  fluctuations  comes  into  play.  Their 
effects  can  be  studied  based  on  the  symmetry-breaking 
pattern  of  the  ordered  phase.  For  a  2D  system  in  the 
continuum,  the  homogeneous  phase  at  d  <  dc  has  con¬ 
tinuous  translational  and  two-fold  rotational  symmetries 
(rotation  by  n)  at  6  ^  0.  For  d  >  dc,  the  stripe  phase 
spontaneously  breaks  the  continuous  translational  sym¬ 
metry  in  one  direction  and  leads  to  one  gapless  Goldstone 
mode  as  required  by  symmetry.  In  a  2D  system  with 
d  >  dc,  a  we  increase  the  temperature  from  T  —  0,  such 
a  Goldstone  mode  will  destroy  the  long-range  order  of 
the  stripe  phase.  However,  a  quasi-long-range  power-law 
correlation  would  still  remain  for  temperature  below  a 
transition  temperature  Tkt  For  T  <  Tkt,  although 
the  real-space  density  oscillation  of  the  stripes  is  hard 
to  observe  due  to  the  lack  of  a  true  long-range  order, 
the  density  profile  in  Fourier  space  will  have  two  peaks 
located  at  ±Q  which  decays  as  a  power-law  function  of 
\qTQ\-  Above  Tkt,  the  correlation  becomes  short-ranged 


and  the  phase  transition  at  T  =  Tkt  belongs  to  the  usual 
Kosterlitz— Thouless  (KT)  university  class. 

For  dipolar  molecules  in  a  2D  optical  lattice,  if  the 
wavevector  of  the  stripe  phase  is  commensurate  with  the 
lattice  wavevector,  the  stripe  phase  breaks  no  continu¬ 
ous  symmetry  and  hence  no  gapless  Goldstone  mode  is 
present.  Therefore,  the  long-range  order  of  the  stripe 
phase  remains  at  finite  T  below  a  critical  temperature 
Tc-  However,  if  the  stripe  phase  is  incommensurate  with 
the  underlying  lattice,  the  KT  physics  described  above  is 
expected. 

In  summary,  we  find  that  the  2D  dipolar  Fermi  gas 
undergoes  a  stripe  ordering  at  large  dipolar  interaction 
strength,  i.e.,  the  density  wave  instability  at  finite  mo¬ 
mentum,  instead  of  phase  separation.  For  the  single 
component  case,  this  instability  occurs  regardless  of  the 
dipole  orientation.  This  arises  from  the  fact  that  the 
Fourier  transform  of  the  dipolar  interaction  increases  lin¬ 
early  with  wavevector.  For  the  multicomponent  case,  due 
to  the  short-range  Hartree  interaction,  the  stripe  insta¬ 
bility  only  exists  at  0  >  Oq. 
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